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ABSTRACT 


A moving -horizon inferential state-estimation technique is 
described which uses simulated 'experimental' data (obtained by 
adding random fluctuations to model results) on temperature and 
viscosity to study bulk polymerization of free radical systems. 
The short-term predictive capability of this technique is found 
to be quite good (except for the case when the curve-fitting 
horizon lies just before the gel-effect region, and does not 
include even a single 'experimental' point from this region). A 
considerable amount of ringing (oscillations between the lower 
and upper bounds) is observed in the values of the estimated 
parameters. The ringing can be reduced significantly by 
narrowing down the range of parameter values , or by taking 
'experimental' data which are spaced further apart. The 
long-term predictions of the model using the estimated parameters 
may or may not be excellent, depending on the 'local' values of 
the parameters. Because of this, periodic use of state-variable 
estimation followed by the determination of the optimal 
temperature history in the future, could be a feasible strategy 
for experimental on-line optimizing control of bulk free-radical 
polymerizations which exhibit significant amounts of the 


Trommsdorff effect. 
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CHAPTER 1 


INTRODUCTION 

In free-radical polymerizations, the termination, 

propagation and initiation reactions (see Table 1) become 

diffusion-controlled as the conversion of monomer increases and 

the viscosity of the reaction mass goes up. The manifestations of 

these diffusion-controlled reactions (sharp and large increases | 

in monomer conversion, x^, with time after a certain period, 

increase in the weight average chain length, jn ; see Nomenclature 

w 

for symbols and definitions) are commonly referred to as the gel j 

1 

1 2 

(or Trommsdorff ) , ' glass and cage effects, respectively. A good ^ 

model is required to account for these effects mathematically, so I 

I 

that it can be used for design, optimization and control j 

[ 

purposes. A sample polymerization system exhibiting these I 
phenomena is that of polymethyl methacrylate (PMMA) , an important | 
commodity plastic. | 

A vast amount of research has been reported in the open ! 

i 

literature in the last two decades on the development of 1 

I 

3 ^ 

theoretical models for MMA polymerization. Chiu et al . have | 

I 

developed a model having a molecular basis and used the 
Fujita-Doolittle free volume theory to account for the 
diffusional limitations of the termination and propagation rate 
constants, k and k . In this model, a set of algebraic 

C It 

equations were written for the gel and glass effects, using the j 
initial nutnber average chain length, as a paraTneter. Later, 

Achilias and Kiparissides^ ' ^ developed a model using the free 


TABLE 1 


KINETIC SCHEME FOR POLYMERIZATION OF MMA 
(BULK AND SOLUTION POLYMERIZATIONS) 


1. Initiation 


2. Propagation 

3 . Termination 

by combination 

by disproportionation 


2R 


k. 

R + M > P 


Pn. M — ^ P^^^ 


P + P — ^ D 
n m n+m 


P + P > D + D 

n m n m 


4 . Chain transfer to monomer 


5. Chain transfer to monomer 
via solvent 


P + M 
n 



P, + D 
1 n 


P + S > S' + D 

n n 

„ . fast „ . _ 

S + M > S + 

n 

or 


P + 
n 


M 


D + P. 


n 


volume theory of Vrentas and Duda^ and accounted for the 
diffusional effects on k , k , as well as the initiator 

L p 

efficiency, f. This model used the initial concentration of the 

initiator, as a parameter. Normally, industrial reactors 

are operated under nonisothermal or semibatch conditions with 

addition and removal of components like initiator, solvent and 

monomer. The models of Chiu et al . and Achilias and 

Kiparissides cannot be applied to these situations since the 

value of ^ or [I]^ is not precisely defined for such cases, 

7 

Ray et al . developed a model which did not have such 
limitations. They assumed the initiator efficiency to be constant 
and obtained expressions for two parameters of the model, ^^(T) 
and 0 (T) , by curve-fitting experimental data®'^ taken in small 
ampoules under isothermal conditions. The predictions made by 
their model for idealized, nonisothermal conditions (step changes 
in temperature) as well as for idealized semibatch operating 
conditions (intermediate addition of a solution of initiator in 
monomer) , were found to be in good agreement with experimental 
observations ' . This confirmed the adequacy of the model for 

simulation of polymerizations under industrially relevant 
conditions of operation. 

12 

More recently, Seth and Gupta modified the model of Ray 

et al."^ and also considered the variation of the initiator 

efficiency, f, with time. Their model used three parameters, 

0 (T) , 0 (T) and 0^ (T) , representing the gel, glass and cage 

t p 

effects, respectively. They again curve -fitted the experimental 
data®'^ on MMA polymerization under isothermal conditions in 
small ampoules (both for bulk polymerization with AIBN 


3 


initiator as well as for solution polymerization in benzene with 

9 

benzoyl peroxide ) . Predictions of their 'tuned' model were found 

to be in better agreement with experimental data^^'^^ taken under 

idealized nonisothermal and semibatch conditions, than the model 

of Ray et al . The complete set of equations for this model and 

the values of several parameters and physical properties 

required, are given in Appendices I and II. ' 

With an appropriate model available for free-radical 

polymerizations, one can now think of studying the on-line 

optimizing control of these processes. An important step in such 

studies is the measurement of the state of the system at 

different times, t. Densitometers and gel permeation 

chromatography have been used in some experimental control 
. 13-15 

studies of solution polymerizations for estimating monomer 

conversions, / and the number average chain lengths, 

[or the weight average chain lengths, u (t) ] . However, these 

experimental techniques cannot be used conveniently for bulk 

polymerizations, and inferential state estimation techniques have 

12 16 

to be resorted to. Seth and Gupta and Chakravarthy et al. 
have suggested the use of experimental values of the viscosity, 
7}(t), of the reaction mass [along with the measured values of the 
temperature, T(t)] for such purposes. Recently, Embirucu et 
al.^"^, have surveyed the open literature on advanced control 
strategies for polymerization reactors, and have found that very 
few studies have been reported on property estimation techniques 
for the control of bulk polymerizations. There is, thus, a 
j^ggd to explore whether model— based inferential state 
estimation using T(t) and T](t) is, indeed, feasible. If so, it 


4 


could be used experimentally (in the future) for on-line 
optimizing control of bulk polymerization reactors. This study 
presents some theoretical work along these lines for the sample 
system, PMMA. 


5 



CHAPTER 2 


FORMULATION 


Model 

The kinetic scheme for MMA polymerization is shown in 
Table 1. A set of ordinary differential equations (ODEs) 
representing the mass balances and moment equations for a 
semibatch reactor can easily be written. These are of the form 

dx / dt = F ( X ) ; X (t=0) = x^ (1) 

where x(t) is the state variable vector defined by 


X = [ I, M, R, S, 

^0 ^1' 

^2/ ^Q, 

“2- Cm' 

^ml 

] (2) 

The k^^ (k = 0, 1, 2) 

moment 

of the 

radical 

and 

dead 


macromolecular species, and respectively, are represented 

by and Hy, (see Nomenclature) . The exact equations are shown in 
Appendix I. Eqs . 1-12 of this Appendix account for intermediate 

additions and removals of components like initiator, solvent and 

monomer through the terms ^Is' ^Im' ^vs ^vm' 

conversion of monomer is defined as 

X = 1 - (M / t 1 ) (3) 

m ml 

where C , is the net monomer added to the reactor since the 
^ml 

beginning of the operation. 

The variation of the rate constants, k^ (= k^^ for PMMA, 

since k^^ = 0), k , as well as the initiator efficiency, f, are 

^ . .12 . 

described by the algebraic Eqs. 13-27 in Appendix I . It is 


6 



observed, that and f depend on the current values of T, 

^n' ^o' ^1' iiot on the initial conditions as in the 

models of Chiu et al . and Achilias and Kiparissides ^ ' ^ . This is 
why the model of Seth and Gupta^^ can be applied easily to 
nonisothermal and semibatch operations of reactors. The values of 
several properties and parameters required to integrate the 
equations in Appendix I are given in Appendix II (for hulK 
polymerization of MMA using AIBN initiator) . The three 
parameters, 6^. (T) , ©plT) and 8^ (T) , have been expressed^^ in 
terms of second order polynomials in 1/T as shown in Appendix II. 
The equations in Appendix I can be integrated (for a given set of 
initial conditions ) using the NAG library program, D02EJF 
(which uses Gear's technique ), to obtain histories of monomer 
conversion and the number and weight average chain lengths . The 

7 

value of the parameter, TOL, required in this code was 10 (the 

results were insensitive to decreases in the value of TOL) . The 

values of the coefficients, a^-a^, and c^-c^, in the model 

12 

have been obtained earlier by curve-fitting experimental data 

under isothermal conditions on bulk® and solution® 

polymerization at several temperatures. In the present work, the 

parameters, a^-a^, re- tuned using only the 

8 

bulk polymerization data since these were of interest . These 
parameters were found to be quite close to values obtained by 
Seth and Gupta^^, and the new values are given in Appendix II. 

An equation relating the viscosity, ij, of the reaction mass 
and the temperature, T, to the state of the system (monomer 
conversion, x^, and weight average molecular weight, M^) is 
required so that model -based inferential state estimation can be 



19 

done, Moritz has suggested the the use of the Lyons- Tobol sky^^ 
equation. This equation is a one-term adaptation of the more 
general series expression relating the ^ specific viscosity^ 
(t//i]soi) "1/ to the product of the intrinsic viscosity, [tq] , and 
the polymer concentration, ^polym* modified the Lyons -Tobol sky 
equation slightly to incorporate a higher-order (quadratic) term 
as shown in Appendix I (Eq. 28) . This was done to obtain better 
predictions for some preliminary experimental data on nonreacting 
PMMA-MMA solutions at 50*^0 and 70*^C, taken on a Haake^ Rotovisco 
RV20 viscometer in our laboratory^^ (see Appendix III) . Such an 
empirical modification also helped in suppressing the temperature 
variation of the parameter, b, in the adapted Lyons -Tobol sky 

equation. The parameter, in this equation, which is a 

1 

function of temperature, was tuned using the (preliminary) 
viscosity data generated in our laboratory, and two coefficients, 
d. and d_, relating k„ to T using a linear variation, were so 
obtained. The values of d^ and d^ are given in Appendix II. It 
must be emphasized that more , experimental data on the viscosity 
of (nonreacting) PMMA-MMA solutions at several additional 
temperatures and concentrations must be taken before the values 
of d^ and (given in Appendix II) can be used for experimental 
studies. However, since the main objective of the present study 
was to explore whether viscosity measurements can be successfully 
used for model-based, on-line inferential state estimation, the 
use of the order-of -magnitude estimates of d^ and d^ given in 
Appendix II (based on the results taken over a restricted range 
of experimental conditions), is appropriate. 



Generation of ^Experimental" Data 

In order to explore the idea of using measured values of 

T(t) and T 7 (t) for inferential state estimation of bulk PMMA 

reactors, one needs to have experimental viscosity data for a 

given T(t) for a polymerizing system. Since this is not yet 

available, (pseudo) 'experimental' data are generated using the 

model itself, and some noise is superposed on the 'smooth' 

model -predictions to simulate actual experimental data. One of 

the optimal temperature histories for bulk polymerization of MMA, 

as provided by Chakravarthy et al.^^ recently, is selected 

(desired value of ii = 1800, and of monomer conversion, x = 

0.94, to be obtained in the minimum reaction time) for this work. 

This 'smooth' temperature history, T (t) , is curve-fitted using 

sm 

t h 

a 17 order Chebyshev series (using the NAG library program, 
E02ADF) so that it can be provided as an input to the 
reactor-simulation program. A noise is superposed on ^^^(t) , 
using a random number, [R(t)], generator to give what can be 
considered as the 'experimental' temperature history, T(t) : 


T (t . ) = T (t . ) 
exp 1 sm 1 


+ [4R(t^) 


In Eq. 4, R(t) is a random number lying between 0 and 1, 
generated using the NAG subroutine G05CCF. The term, 4R(tj|^) - 2, 

is used so as to give the amplitude of the noise as 2°C. The 
'experimental' data points are generated using Eq.4 with = 

i(At), where At is taken as 0.5 min. The viscosity of the 
reaction mass for bulk polymerization of MMA has then been 
generated using the initial value of the initiator concentration, 
[I]^, as 15.48 mol/m^, and using Tg^p(tj^} [with linear 


9 



interpolation used between consecutive points of T (t.) in 

exp 1 

D02EJF] . It is observed that the model values of T)(t) so 
obtained are quite smooth and we must, therefore, add on a noise 
aqain so as to provide simulated 'experimental' viscosity data. 

The following equation is used for this purpose : 

^exp^^i^ " R(t^)] T){t^) ; t^ = i(At) (5) 

Use of Eq . 5 leads to fluctuations in ranging from O.Stj to 

1.5t], a fairly large range as compared to the possible errors and 
fluctuations in real experimental viscosity data. 

The 'experimental' temperature and viscosity data so 

' ' I 

generated (by simulation) are shown in Figs. 1 and 2. These data 
are used for studying model-based inferential state estimation. 

The final values of x^ and using T(t) given by the solid curve 
in Fig. 1 are 0.9408 and 1748, respectively (these differ 
slightly from the values in Ref. 16 because our parameters are 
slightly different) . 

Inferential State Estimation 

In a typical polymerization reactor in which on-line i 
optimizing control is implemented, one would have available, at j 
time t = t^, a set of experimental values of temperature and 

viscosity, [i(At)] and /nexp ; i = 0, 1, 2,..., 

t^/(At). A short time-horizon for curve -fitting, t^- N^(At) t 
5 t^, is selected and the experimental points lying (only) in 
this horizon are used to estimate the values of several of the 
parameters in the model. Sequential quadratic programming 
(SQP)^^'^^ has been used in this work to obtain best-fit values 



of four parameters, a^, b^, and (referred to as a^ - d 3 _; 
while assuming all the remaining parameters to be the same as 
given in Appendix II — these being called reference values 
henceforth) . The objective function used for the optimization is 


Min E 






^exp^tj ^ - ^th^^i) - 


N, 


^th^4^ 


( 6 ) 


where is the experimental value of the viscosity, and is 
the value predicted by the model corresponding to the values of 
^ 2 ," order to obtain the model values of the several 
state variables (see Eq. 2) and for t^- N^(At) s t :£ t^ (as 
required in Eq. '6) , the initial conditions are taken to be the 
model values at t = t^- N^(At) as' computed in the previous 
iter:ation (and stored). Also, a constant temperature is assumed 
for the duration t^- (At) ^ t t^, which is the average of the 
'experimental' values of temperature during this period. The NAG 
subroutine, E04UPF using SQP, has been used for obtaining the 
best -fit values of a^- d^ (which minimize the sum of square 
error, E, in Eq. 6) . Table 2 gives the values of the 
computational . parameters used with this code. The SQP procedure 
also requires initial guesses and bounds for the parameters, a^- 
d^, for each value of t^. The initial guesses supplied for the 
first iteration [t^= N^(At)] are .given in Table 2. For 
subsequent values of t^, the optimal values of the parameters in 
the previous iteration are used as the initial guess. The bounds 
on the parameters are also listed in Table 2. The bounds on ~ 
c^^ have been chosen .such that variations of al^ut a decade (at 
any given temperature) are peimaitted in the values of 0^ and 


TABLE 2 


INITIAL GUESSES AND BOUNDS FOR GENERATING THE SOLUTION IN THE 


FIRST 

CURVE-FIT HORIZON, AND 

PARAMETERS USED IN 

THE SQP CODE 

Parameter Initial Guess 

Lower Bound 

Upper Bound 


1.2416x10^ 

1.2346x10^ 

1.2446x10^ 


8 . 0673x10^ 

8 . 0071x10^ 

8 . 1071x10°’ 


2.0168x10^ 

2.014x10^ 

2 . 024x10^ 


0.3118 

0.309x10° 

0.315x10° 

^2 

1.0314x10^ 

1.03018x10^ 

1.0334x10^ 

^2 

7 . SXIO'^ 

7.4845x10^ 

7.5168x10^ 

^2 

1.445x10^ 

1.4524x10^ 

1.4556x10^ 

^2 

9 . 93x10"'^ 

9 . OxlO'^ 

1. OxlO"^ 

^3 

2.2735x10^ 

2.2669x10^ 

2.27743x10°^ 

^3 

1. 765x10"^ 

2 . 75955x10"^ 

1.76999x10°^ 

^3 

2 . 7x10"^ 

2 . 6979xlo'^ 

2 . 7083xlo'^ 


Parameters 

in SQP Code^^ 



Parameter 

Value 



10 


10 

func^^ 

O 

H 

1 

O 

H 

pr 

^^^•act 

10 


(note that a^- are independent of temperature) . The bounds 
on d^ are chosen so as to permit a fairly large variation in the 
viscosity of the reaction mass. Much smaller bounds would 
normally be used in experimental on-line control work. 

In addition to obtaining best-fit values of - d^ in any 
iteration of curve-fitting (each such iteration requiring several 
iterations in the SQP code for convergence to be attained) , we 
have used the reactor simulation code to predict values of the 
state variables and of T?,-j^(t) for 'future' times, t^ :£ t :£ t^ + 
N 2 (At) [using '^expt^^^ this prediction horizon]. These 
theoretical predictions could be compared to the 'experimental' 
data points of viscosity (in the prediction horizon) to find out 
how good is the (short-term) predictive capability of the 
model -based state estimation technique described herein. The 
computer code starts from t = 0 and uses the first points 
(first curve-fitting horizon) for curve-fitting (followed by 
prediction of N 2 points in the future) . This completes the first 
iteration of curve -fit ting . The experimental values of viscosity 
for 0 + (At) < t :£ (N^ + N^) (At) are then used for the second 
iteration (second curve-fitting horizon) . This continues till 
the end of polymerization. The values of can be selected so 
as to have some overlap between the successive curve-fitting 
horizons. 

A few checks were made on the computer code to ensure that 
it is free of errors. When the optimization program was run 
using all the 'experimental' points (with no fluctuations 
introduced) and initial guesses were provided which differed 
from the reference values (Appendix II) of a^- d^. 


the 



optimization code converged to the reference values. Some cnher 
checks were also conducted and led to the conclusion that the 
code was relatively free of errors. The CPU time taken for the 
optimization problem was about 3 0 s on a suner-mini HP 9000/73 5 
mainframe computer. 



CHAPTER 3 


RESULTS AMD DISCUSSION 

The optimization program was run on the 'experimental' data 
of viscosity and temperature (Figs. 1 and 2), using the following 
parameters (called reference values again) 

= 5 
= 5 
= 1 

= 0.5 min ( 7 ) 
Thus, five 'data' points at a time are taken for curve-fitting, 
and the curve- fitting horizon shifts by one 'experimental' point 
at a time in this reference run (with four points common between 
consecutive curve-fitting horizons) . 

Fig. 3 shows the piece-wise fitted curve of t] {t) alongwith 
the 'experimental' data points. The solid curve in Fig. 3 
consists of segments which extend for one At period (the 
non-overlapping domain between two consecutive curve-fitting 
horizons) . The fit of the viscosity 'data' is seen to be quite 
good. Fig. 4 shows the variations of the parameters, - d^, 
from horizon to horizon (the end-values of t of the curve-fitting 
horizons are shown in Fig. 4 on the abscissa) . It is observed 
that there is a considerable amount of 'ringing' in the values 
of these parameters between the lower and upper bounds. It is 
demonstrated later in this paper that this is because we are 
trying to curve- fit only a. small number (5) of experimental data 


N. 


N, 


^3 

At 










reference case. 'Experimental' viscosity data, shown by 





case, corresponding to the curve-fit of Fig. 






points in any one iteration, and that the ringing can be reduced 

by changing the bounds of the four parameters, as well as by 

increasing At. The corresponding variations in 6 (T) , 0 (T) , 

^ P 

0 ^ (T) and are shown in Fig. 5. The ringing in the values of 

~ reflected as similar sharp fluctuations in the values 

of 0^, 0p, 0^ and kjj. In addition, there is a general change in 

these parameters associated with their temperature dependence^^ . 

The monomer conversion, M and M for the reference case are 

n w 

shown in Figs. 6 and 7. It is observed that the model 

predictions agree quite well with the 'experimental' points 

corresponding to the temperature 'data' of Fig. 1 (it may be 

noted that the viscosity 'data' were generated by introducing 

additional fluctuations to values computed using the 'data' of 

Figs. 6 and 7) . The final values (at t = 30 min) of x and 4 

m n 

are found to be 0.9307 and 1743 respectively, for the 

curve-fitted case, as compared to the values of 0.9408 and 1748 
corresponding to the temperature history shown by the solid 
(smooth) curve in Fig. 1. The small disagreements in these 
values are because of fluctuations introduced in T(t) and T?(t) . 

The optimal parameters obtained in any curve- fitting horizon 
are used with the temperatures in the corresponding 

prediction-horizon (five further points, since = 5) to predict 
the viscosity. The short-term prediction capabilities are shown 
in Fig. 8 . It is clear from this figure that if the 
curve- fitting horizon is far ahead of the gel-effect, the 
short-term (N 2 = 5) predictions are quite good. When the 

curve - fitting horizon lies just before the gel-effect, the 
short-term future predictions need not excellent. The 
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CD _ tn 
cs 
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(T) , corresponding to the 













predictions, however, improve considerably when even a single 

data point from the gel-effect region is included in the 

curve-fitting horizon. This is shown in Fig. 8b. Similar 

results are obtained for short-term predictions of x ii and tz . 

m n ^w 

All these reference results indicate that 'experimental' data on 
T (t) and 7) (t) can indeed be used with a model to estimate the 
state variables of the system, e.g., x^(t) , , h^(t) and 

7] (t) . 

It may be added that long-term predictions using the present 
technique may not be as good. For example, use of five 
'experimental' points of 7}(t) for curve-fitting the values of a^ 
- d^, and their subsequent use in the model with T(t); t 2: t^, 
did not lead to as good an agreement of the predicted 7 )(t) with 
'experimental' values in the gel-effect region (see Figs. 9a and 

10a) . However, use of another set of five points leads to 

reasonable long-term predictions as shown in Figs. 9b and 10b 
(this is a consequence of the ringing in the values of the 

parameters) . This may pose seme problems in the use of this 

technique for on-line optimizing control purposes with 
minimization of the reaction time as the objective function. We 
feel, however, that multiple use of this technique (with shifting 
curve- fitting and prediction horizons) may overcome this problem, 
and the computed optimal temperature history may be satisfactory. 
Work is in progress to see if, indeed, this conjeccure is true. 

Having established the applicability of the technique 
described herein for inferential state-variable estimation, we 
now study the effect of varying several of the other parameters 
involved. Instead of using the most recent optimes.! values of 








case) for two curve-fitting horizons. Notation 






as the initial guess values for any curve-fitting horizon, 
we used the initial values given in Table 2 for all the horizons. 
There was no perceptible change in the fit of the 7 ?(t) data. 
Almost the same degree of ringing of the parameter values was 
observed. The CPU time was also about the same (increased by 
about Is) . These and other results not included in this paper 
for reasons of brevity, can be supplied on request. 

In the reference run, we assumed a^, b^/ b^ , c^, c^ and 
*^2 constants (at the values given in Table 2) , and we 
obtained best -fit values of a^, b^, c^ and d^ . We explored if 
there would be any improvement if all the 11 parameters (a^ - d^) 
were used for curve-fitting the 7?(t) 'data'. The bounds of the 
remaining 7 parameters provided to the SQP code, as well as the 
initial guesses, are also given in Table 2. The bounds on a^, 
a^ / b 2 , b^ , and were selected so that each of them 
individually could give variations of about a decade in the 
values of 0^, 0^, and 0^. The CPU time increased by about four 
fold, to 111 s. We observed that there was not much improvement 
in the fit of the viscosity data. Also, the degree of ringing in 


the 11 parameter, did not decrease much. Much larger 
oscillations were observed in the plots of 0^ (T) , 0p(T) and 0^ (T) 
for this case, possibly because of the adding-up of the effects 
of each of the additional parameters. Use of 11 parameters is, 
thus, not recommended since it provides no additional advantages 
over the use of 4 parameters (reference case) . 

We also studied the effect of narrowing the bounds on a^ - 
d^ by about 50%. The final fit of the viscosity data was about 
as good as that observed for the reference run (Fig. 3) . 


But 



as shown in 


the ringing of the parameters, - d^, reduced 

Fig. 11. No worsening was observed in the short-term prediction 

capabilities. The final values of x , M and M were observed to 

m' n w 

be close to the reference values. The CPU time came down to 10 
s. It can, thus, be inferred that narrowing the hounds for a^ - 
d^ is desirable for experimental, on-line optimizing control, and 
one needs to make a few test-runs to obtain the best-values. 

We next studied the effect of increasing the value of N^, 
the number of points used in the curve -fitting horizon, from 5 to 
8 . The fit of viscosity was observed to be as good as that shown 
in Fig. 3. The fluctuations in - d^ reduced slightly as 
compared to the reference case, as shown in Fig. 12. The final 
values of x^, and were in the desired range, as in the 

reference case. The prediction capabilities were also found to 
be about the same as for the reference case. The CPU time came 
down, as expected, to 15 s. 

We, then, studied the effect of increasing the time 

interval. At, between 'data' points to 1 min (reference value = 
0.5 min). The alternate 'data' points from the reference case 
(Figs. 1 and 2) were used for this purpose. The results obtained 
showed equally good fit of the viscosity data, as for the 

reference' case. Fig. 13 shows that the degree of ringing in a^ - 

d^, with At = 1 min, is much less. The prediction capabilities 

and the final values of x^, and M^, were found to be as good 

as for the reference case. The CPU time was 10 s. It is, thus, 
observed that larger values of At and narrower bounds on a^ " 
than used for the reference case, are to be used for on-line 

The best values of these can be 


optimizing control work. 












obtained only after combining the present code with a program 
using genetic algorithm^^ to obtain optimal future temperature 
histories periodically. 

We then explored the reasons for getting ringing in the 
values of the parameters. Isothermal (50°C) 'data' on 7 ?(t), 
without any fluctuations introduced, was taken. This is shown in 
Fig. 14. The entire set of 66 points were curve-fitted to yield 
single values of the four parameters, a^ - d^, as shown bycurves 
A in Fig. 15. The same 'data' was curve-fitted using 35 data 
points each, with 4 overlapping points (see Table 3) . Curves B in 
Fig . 15 show sudden and significant changes in the value of at 
least one parameter, c^. Cases C - B in Fig. 15 show a 
significant amount of ringing as the number of curve-fitting 
horizons is increased. The fit of viscosity is equally good (as 
shown in Fig. 14 by the solid curve) in all the cases. It is, 
thus, clear that ringing is introduced by the use of only a few 
data points of viscosity, at a time, and that it is not caused by 
fluctuations in the data. The long-term predictions (for = 5), 
again, could be excellent or poor, because of the ringing in the 
values of the parameters . 




Fig. 14 Experimental ' viscosity data (points) under isothermal 
(50 °C) conditions with (I] = 15.48 mol/m^, in 












TABLE 3 


DETAILS OF PARAMETERS USED FOR FITTING ISOTHERMAL VISCOSITY 

OF FIG. 14 (At = 5 min) 


Curve 


^3 

A 

66 

- 

B 

35 

01 

C 

40 

02 

D 

05 

01 


DATA 



CHAPTER 4 


CONCLUSIONS 


A model-based, inferential state-variable estimation 

technique has been described which uses a few of the most recent 

'data' points on the viscosity of the reaction mass (moving 

curve- fitting horizon) to predict x , and u . Short-term 

m n w 

predictions in the moving prediction-horizon are found to be 
quite good, except when the curve-fitting horizon is very near to 
the get -effect region and does not include any data point in this 
region. Long-term predictions may or may not be excellent, 
depending on which curve-fitting horizon is used. Multiple 
curve- fitting followed by periodic determination of the optimal 
temperature history in the remaining period of polymerization 
could provide a satisfactory scheme for on-line optimizing 
control . 
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CHAPTER 5 


SUGGESTIONS FOR FUTURE WORK 

Combining this work, periodic state-variable estimation, 
with the determination optimal temperature history in the future 
it could be seen whether this could be a feasible strategy for 
experimental on-line optimizing control of bulk free-radical 
polymerizations exhibiting Trommsdorff effect. 
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APPENDIX I 


MODEL EQUATIONS FOR MMA POLYMERIZATION IN SEMIBATCH REACTORS^^ 
(BULK AND SOLUTION POLYMERIZATION*) 
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^td " 2 '"tc' V, 
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Appendix I (contd. . . .b) 


dju, 


dt 




dll^ X XX 'x^ 

SE^ ■ <V + + k,— °„ + k J- 

ac s t t ^ ^tc 
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11, 


m 
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^"Ira'- ■ -vm 
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Appendix I (contd....c) 
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= k° exp (-E^/RT) 
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28 . T) 


■^sol_^ Solym (■ 
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30. E-n] = K m; 
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* no inert solvent (benzene) used in this work 

Note: variables not defined in the Nomenclature are defined 


Refs. 7 and 12. 
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APPENDIX II 


PARAMETERS USED FOR BULK POLYMERIZATION OF MMA WITH AIBN 

“ 966.5 - 1.1 (T - 273.1) kg/m^ 

= 1200 kg/m^ 

f = 0.58 

0 

k° = 1.053x10^^ s"^ 

a 

k° » 4.917x10^ m^/mol-s 

P t ^ 

k°j * 9.8x10^ m^/mol-s 

td,o 

k^ = 0 . 0 

tc 

kg - 0.0 

k. » k„ 

1 p 

k - 0.0 

8 

» 128.45 kJ/mol 

E = 18.22 kJ/mol 

P 

Eg^ * 2.937 kJ/mol 

(MW )» 0.10013 kg/mol 

tn 

(MWj)« 0.06800 kg/mol 

Parameters for the Cage, Gel and Glass Effects 

9.13x10“^ m^/kg 
8.22x10”^ ra^/kg 


Contd . . 


Appendix II (contd. . . b) 

V* = 7.70x10“^ m^/kg 

M. = 0.18781 kg/raol 

J IT 

r = 1 

= 0.149 + 2.9x10'^ T(K) - 273.lj 

Vfp = 0.0194 + 1.3x10"^ T(K) - 273.1 - 105 ; 

for T< (105+273 .1) K 

Mark-Houwink Constants for Intrinsic Viscosity (see Eq. 29, App. I) 

K = 6.75 X 10'^ m^/kg 

a = 0.72 

K and a assumed to be (almost) independent of T 

Parameters for the Adapted Lyons -Tobolsky Equation (Eq. 28, App. I) 
+ d^T 

d^ = 0.3118; d 2 = 9 - 93 X 10 K ^ (this work) 

- 3 3 

b = -3.5x10 m /kg (assumed independent of T) 

-I C QO Q 

^sol "" 0.099 + 496/T) / t Pa-S (Ref. 5) 

Correlations Used For Curve - Fitting 

logio [ , s] = a^ - a^d/T) + a 3 (l/T^) 

log^O^^p^T)' = ^1 ■ b^d/T) + b3(l/T^) 

log^Q [1O^0£ (T) , m^mol'^] = - c^d/T) + c^(l/T^} 

Contd . . . . c 
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Appendix II (contd. . . c) 



= 1.2408x10^ 

; 

= 1.0314x10^ 

; 

= 2. 2735x10"^ 



= 8.0593x10^ 

'• ^2 

= 7.5xl0'^ 

'• ^3 

= 1.765xlo'^ 



= 2.016x10^ 

' ^2 

= 1.455x10^ 

; 

= 2.70xl0’^ 

(this work) 
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VISCOSITY OF NONREACTIMG MMA - PMMA SYSTEMS 


No. 

Cone . of polymer 
% by weight 

T,°C 

Ti, Poise 

1 

38.87 

50 

2800 

2 

38 . 87 

70 

2500 

3 

31.60 

50 

200 

4 

31.60 

70 

120 

5 

20.23 

50 

6.5 

6 

20.23 

70 

2.6 

Since 

M were not measured, 

w 

these 

were obtained by 


computation, as values corresponding to a reacting system at the 

3 

same T and polymer concentration; C13 q = 15.48 mol/m . 


